Photoluminescence of nanocrystals embedded in oxide matrices 
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We used the theory of finite periodic systems to explain the photoluminescence spectra dependence 
on the average diameter of nanocrystals embedded in oxide matrices. Because of the broad matrix 
band gap, the photoluminescence response is basically determined by isolated nanocrystals and 
sequences of a few of them. With this model we were able to reproduce the shape and displacement 
of the experimentally observed photoluminescence spectra. 

I. INTRODUCTION 

ZnO-Gc and Si02-Ge nanocomposites are intensively studied for their possible use as optoelectronic devicesii2i2i2£. 
A common method for preparing these systems is by alternate deposition of oxides and Ge, followed by an annealing 
process to increase the diffusion of nanocrystalites (NC) into the adjacent layers^iSiS. This process leads also to 
homogenize the sizes and spacial dispersion of NCs. It is widely known that the photoluminescense-spectra shift is 
strongly related to NC's size. Although there have been several attempts to understand the relation between the 
morphology and the photoluminescence (PL) response of embedded NC in oxide matricesA^, none of them lead to a 
simple and precise explanation. The purpose of this work is to understand and to explain this relation from a simple 
model and a few physical parameters. 

To calculate the optical response of these rather complex systems, we use both the bona fide multiple-quantum-well 
electron and hole eigenvalues and eigenfunctions, obtained from the theory of finite periodical systems^ (TFPS) and 
a simple assumption, experimentally supported, on the statistical distributions of NC sizes. We show that the main 
features of the PL spectra (in Ref A for Ge doped ZnO films annealed at several temperatures) , can clearly be identified 
in terms of the interband excitonic transitions calculated for multiple quantum well (MQW) systems, whose well and 
barrier parameters are drawn from the experimental system. The agreement with other kind of experimental results, 
as those reported by by S. Takeoka et alA for Ge NC in S1O2, is also good when the nanocrystals-size dispersion is 
taken into account. 



II. THE MODEL 
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In Fig. ^ we observe micrographs of Ge cristalites embedded in a ZnO matrix. At low NC densities we have 
almost isolated NCs (see Fig^ (a)), while at higher NC densities one can find groups of NCs aligned and almost 
equidistant (see Fig. H(b)). Although a single NC can be better modeled as a quantum dot, the relevance of the 
longitudinal coherence in the energy spectra of a sequence of NCs, leads us to model the sequence as a multiple 
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FIG. 1: Ge cristalites embedded in a ZnO matrix for (a) low NCs densities and (b) high densities. 
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FIG. 2: (a) Finite periodic square-well modulation of the valence and conduction band edges along a set of n QW's. Also we 
show the wave functions ^1 x and tyf 1 for the first energy level of the first subband in the conduction and valence band. For 
these wave functions the electron and hole are localized inside the QW's. (b) Room temperature PL spectra of Ge doped ZnO 
films annealed at several temperatures and theoretical peaks position (black arrows). We also show the corresponding NC's 
average diameter used for the theoretical model. 



ID-quantum well potential. Each of these sequences is then treated as a finite periodic square-well potential (see Fig. 
0(a)), and the whole nanocomposite PL response results from the contribution of an ensemble of MQWs, which differ 
basically in the number of NCs, well-widths and barrier-thicknesses. 

In Fig. |3(a) we have a MQW potential for a given set of parameters. We show the energy spectra (in the conduction 
and valence bands), and the square of the wave function amplitude ^1 1 f° r the first energy level in the first conduction 
(valence) subband. In the TFPS, these quantities can easily be calculated. The eigenvalues are obtained fromiS 

(a n e lka -a* n e- lka ) + (3 n -(3* n = 0, (1) 



with k = y2mE /h 2 and m the electron or hole effective mass. The wave functions at any point z in the (j + l)-th 
cell (for a system of length L — nl c ) and energy E are given by 



E) = A[{a p + 7p )K- - f3,%±§-) + (P P + W* - (2) 

' fin ' i n 



aj, Pj are elements of the transfer matrix M(j7 c , 0; E) and a p , f3 p , 7 P , 5 P are elements of a partial cell transfer matrix 
M p (z, jl c ; E). fy L (z, E^ v ) (with l = c,v for conduction and valence bands), is the eigenfunction corresponding to the 
eigenvalue E^ „. To evaluate the PL response we use the following function 

v P^wfjE^) [i-f(E;, ul )} 

X ^ , (w s -E- u + El, v ,y + Ti ■ { > 

f{E[ l „) the Fermi-Dirac carrier distribution functions in the conduction or valence band. 



III. RESULTS AND DISCUSSION. 



In Fig. |2] (b) we present the experimental PL spectra of Ge-doped ZnO films for different annealing temperatures, 
reported in Refi The NC's average diameter increases with the annealing temperature 8 . The peaks position in these 
PL curves can be described in terms of the energy spectra of a single sequence of NCs. For this calculation we use the 
potential and gap parameters corresponding to the T point of ZnO and the annealing temperature, i.e. mg e = 0.12?n e , 
m Ge = 0.23m e , m e Zn0 = 0.24m e , m h Zn0 = 0.45m e , Ec e — 0.66eV, Ez n o — 3.31eV. The black arrows point out the 
theoretical predictions for different possible recombination transitions. In this figure, the assignment [i — > fj, refers to 
the energy transitions E c ^ v — > E^ v . It is well known that increasing the well width a, the subbands move towards 
lower energies while their number increases. We notice in the the experimental PL curves that increasing a the small 
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FIG. 3: Theoretical PL spectra of Ge doped ZnO films for (a) several NCs diameters and (b) barrier widths. 




FIG. 4: At the left the measured PL response for different average diameters of NC-Ge embedded in SiC>2- In the middle our 
theoretical curves, and at the right image we plot the PL peak energy versus average diameter of NC-Ge for: experiment-, the 
theoretical prediction with size-dispersion and without size- dispersion. 

peaks shift toward larger wave lengths, while the main PL peak displaces in the opposite direction. This behavior 
is even much clear in the experimental results of Takeoka eta al. in Ref^, where the secondary-peaks structure is 
absent. 

In these systems the experimental curves arise from the contribution of a large number of NCs which sizes are 
statistically distributed 8 . The PL aspect and peaks relative- height depend not only on the potential parameters, but 
also on their statical distribution. The PL response of a single MQW in Fig. shows the red-shift with a and the 
high PL sensitivity to small barrier-thickness variations. 

To explain the experimental results reported by Takeoka et al. in Hem for Ge Ncs in a SM2 matrix, we need to 
take into account also the NCs size distribution. We assume, according with the experimental results, that the size 
distribution is described by a log-normal function. We then calculate the PL spectra of several sets of NC diameters, 
pondered by their statistical weight. Since the PL peaks position is less sensitive to the number of wells and to the 
NC spacing, we assume for this case a fixed barrier width and also a fixed number of cells. 

At Fig. U|(a) we present the experimental results from Ref. 4 . In Fig. 01(b) we show our calculated curves assuming up 
to 50 contributing MQW sequences with the following parameters: m e Si02 — 0.5m e , m h Si0<2 = 0A9m e ,Esio 2 = 9.3eV, 
corresponding also to the T point of Si02- At Fig. 0] (c), we plot also the peak position versus the NCs average 
diameter. This averaging procedure improves the agreement with the experimental PL spectra. The size-dispersion 
shifts the PL peaks toward lower energies. This is a consequence of the asymmetry of the log- normal distribution. 

In our calculation we neglected e-h Coulomb interaction corrections. In any case, they will produce small red shifts 
in the secondary PL peak positions and will be washed out in the main peak behavior. 



4 



IV. CONCLUSIONS 



To explain the photo-luminescence response of NCs embedded in oxide matrices as function of the average NCs 
diameter, we modeled these systems as a collection of 1-D multiple-quantum-well potentials. Using the finite periodical 
systems theory and a simple NC-size distributions assumption, we shown that the calculated 1-D MQW subband 
transitions, describe correctly the secondary PL peaks behavior, observed by Fan et ali, as function of the NCs 
diameter. We have shown also that the main PL peak blue-shift behavior, observed by Takeoka et al. in Ref. 4 , is 
basically determined by the statistical dispersion of the NCs diameter. 
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